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Abstract Shearlet systems have so far been only considered as a means to ana- 
lyze L^-functions defined on M?, which exhibit curvilinear singularities. However, 
in applications such as image processing or numerical solvers of partial differen- 
tial equations the function to be analyzed or efficiently encoded is typically defined 

PU ' on a non-rectangular shaped bounded domain. Motivated by these applications, in 

(~| , this paper, we first introduce a novel model for cartoon-like images defined on a 

bounded domain. We then prove that compactly supported shearlet frames satis- 
fying some weak decay and smoothness conditions, when orthogonally projected 
onto the bounded domain, do provide (almost) optimally sparse approximations of 
elements belonging to this model class. 



^ ■ 1 Introduction 



It is by now well accepted that L^-functions supported on the unit cube which are 
(^ • C^ except for a C^ discontinuity curve are a suitable model for images which are 

^^ , governed by edges. Of all directional representation systems which provide op- 

timally sparse approximations of this model class, shearlet systems have distin- 
guished themselves by the fact that they are the only system which provides a unified 
.J treatment of the continuum and digital setting, thereby making them particularly 

r> ' useful for both theoretical considerations as well as applications. However, most 

j^ ■ applications concern sparse approximations of functions on bounded domains, for 

instance, a numerical solver of a transport dominated equation could seek a solution 
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on a polygonal shaped area. This calls for shearlet systems which are adapted to 
bounded domains while still providing optimally sparse expansions. 
In this paper, we therefore consider the following questions: 

(I) Which is a suitable model for a function on a bounded domain with curvilinear 

singularities? 
(II) What is the 'correct' definition of a shearlet system for a bounded domain? 
(Ill) Do these shearlet systems provide optimally sparse approximations of the model 
functions introduced in (I)? 

In the sequel we will indeed provide a complete answer to those questions. These 
results push the door open for the usability of shearlet systems in all areas where 2D 
functions on bounded domains require efficient encoding. 



1.1 Optimally Sparse Approximations of Cartoon-Like Images 

The first complete model of cartoon-like images has been introduced in HI, the basic 
idea being that a closed C" curve separates smooth - in the sense of C^ - functions. 
For the precise definition, we let p : [0, 27r] — > M+ be a C^ radius function and define 
the set B by 

B = {x^R^: ||,y||2<p(0),.y= (||,y||2,0) in polar coordinates}, (1) 

where 

sup|p"(0)|<v, p<po<l. (2) 

This allows us to introduce STAF?'{v), a class of sets B with C^ boundaries dB and 
curvature bounded by v, as well as <^^ ( v), a class of cartoon-like images. 

Definition 1 (lH). For v > 0, the set STAR?-{v) is defined to be the set of all B C 
[0, 1]^ such that B is a translate of a set obeying ([T]) and ^. Further, S'^{v) denotes 
the set of functions / on M^ with compact support in [0, 1]^ of the form 

f = f0 + flXB, 



where B e STAR^{v) and /o,/i e C^(E^) with compact support in [0, 1]^ as well as 
E|a|<2l|£>"/dl~<l for each/ = 0,1. 

In H, Donoho proved that the optimal rate which can be achieved under some 
restrictions on the representation system as well as on the selection procedure of the 
approximating coefficients is 

\\f-fN\\l<C-N-^ asA?->- 

where /n is the best A^-term approximation. 
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1.2 Shortcomings of this Cartoon-Like Model Class 

The first shortcoming of this model is the assumption that the discontinuity curve 
is C^. Think, for instance, of an image, which pictures a building. Then the frames 
of the windows separate the dark interior of the windows from the presumably light 
color of the wall, however this frame is far from being C^. Hence, a much more 
natural assumption would be to assume that the discontinuity curve is piecewise C^. 

The second shortcoming consists in the fact that the function is implicitly as- 
sumed to vanish on the boundary of [0,1]^. More precisely, even if the function 
f = fo+ fiXB is non-zero on a section of positive measure of the boundary dB, 
this situation is not particularly treated at all. However, reminding ourselves of the 
very careful boundary treatment in the theory of partial differential equations, this 
situation should be paid close attention. Thus, a very natural approach to a careful 
handling of the boundary in a model for cartoon-like images seems to consist in 
regarding the boundary as a singularity curve itself. 

The third and last shortcoming is the shape of the support [0, 1]^ of this model. 
Typically, in real-world situations the domain of 2D data can be very different from 
being a rectangle, and even a polygonal-shape model might not necessarily be suffi- 
cient. Examples to support this claim can be found, for instance, in fluid dynamics, 
where the flow can be supported on variously shaped domains. In this regard, a suit- 
able model situation seems to be to allow the boundary to consist of any piecewise 
C^ curve. 



1.3 Our Model for Cartoon-Like Images on Bounded Domains 

The model for cartoon-like images on bounded domains, which we now define, will 
indeed take all considerations from the previous subsection into account. For an 
illustration, we refer to Figure [T] 




Fig. 1 Example of a function / belonging to our model class S^ j^{Q.) of cartoon-like images on 
bounded domains. 
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We first introduce STAR-^{v,L), a class of sets B with now piecewise C' bound- 
aries dB and curvature on each piece bounded by v. This will serve us for both 
modeling the bounded domain as well as modeling the discontinuity curve. For this, 
let L e Z+ denote the number of C^ pieces and let v > be an upper estimate for the 
curvature on each piece. Then B £ STAR^{v,L), if B is a bounded subset of [0, 1]^ 
whose boundary dB is a simple closed curve and 

L 

1=1 

where each curve p,- is parameterized by either xi = £',(x2) or X2 ~ Ei{xi ) and £,■ e 

C^{[ai,bi]) such that 

// 
max max LE; < v. 

This allows us to introduce a model class of cartoon-like images in bounded do- 
mains. In accordance with modeling functions on bounded domains, we now con- 
sider functions defined on [0, 1]^; its 'true' domain is brought into play by requiring 
these functions to be supported on 12 C (0, 1)^, which we model as piecewise C^ 
bounded. This ensures that we treat dQ as a singularity curve, which would not 
have been possible when defining the model on Q itself. 

Definition 2. For v > and L e Z+, let f2,B e STAR^{v,L) be such that B C f2°, 
where Q.° denotes the interior of the set ^, and Q. C (0, 1)". Then, (?^^(X2) denotes 
the set of functions / on [0,1]^ with compact support in Q. of the form 

where /o,/i G C^([0, 1]^) with compact support in Q. and L|a|<2 ll^"/iil«> f^ 1 for 
each ! = 0, 1 . 

Later it will become important to analyze the points on boundaries of sets in 
STAR^{v,L), in which the boundary is not C^. For these points, we will employ the 
following customarily used notion. 

Definition 3. For v > and L e Z+, let B e STAR^{v,L). Then a point xq G dB will 
be called a corner point, if dB is not C^ in xq. 

Since the model S^ ^{Q.), while containing the previous model (§y as a special 
case, is considerably more complicated, we would like to make the reader aware of 
the fact that it is now not clear at all whether the optimal approximation rate is still 

||/-/iv||2 <C-A^"^ asA^^oc. 
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1.4 Review of Shearlets 

The directional representation system of shearlets has recently emerged - a first in- 
troduction dates back to 2005 in ITT] - and rapidly gained attention due to the fact 
that, in contrast to other proposed directional representation systems, shearlets pro- 
vide a unified treatment of the continuum and digital world similar to wavelets. We 
refer to, e.g., QEl for the continuum theory, lfT6l l6l [T8l for the digital theory, and 
(SI m for recent applications. Shearlets are scaled according to a parabolic scaling 
law and exhibit directionality by parameterizing slope by shearing, the later being 
the secret which allows the aforementioned unified treatment in contrast to rotation. 
Thus shearlets are associated with three parameters: scale, orientation, and position. 
A precise definition will be given in Section|2] 

A few months ago, the theory of shearlets focussed entirely on band-limited 
generators although precise spatial localization is evidently highly desirable for, 
e.g., edge detection. Recently, motivated by this desideratum, compactly supported 
shearlets were studied by Kittipoom and the two authors. It was shown that a large 
class of compactly supported shearlets generate a frame for l}{M?) with control- 
lable frame bounds alongside with several explicit constructions ifTOl . By the two 
authors it was then proven in Ill51 that a large class of these compactly supported 
shearlet frames does in fact provide (almost) optimally sparse approximations of 
functions in S^ in the sense of 

\\f-fN\\2<C-N-^-{\ogNf asA?^-. 

It should be mentioned that although the optimal rate is not completely achieved, the 
log-factor is typically considered negligible compared to the A'^^-factor, wherefore 
the term 'almost optimal' has been adopted into the language. 



1.5 Surprising Result 

We now aim to analyze the ability of shearlets to sparsely approximate elements 
of the previously introduced model for cartoon-like images on bounded domains, 
S'yi^{Q). For this, we first need to define shearlet systems for functions in L^{Q). 
Assume we are given a (compactly supported) shearlet frame for L^(R^). The most 
crude approach to transform this into a shearlet system defined on L?{Q), where 
Q G STAR^{v,L), is to just truncate each element at the boundary of i2. Since it is 
well known in classical frame theory that the orthogonal projection of a frame onto 
a subspace does not change the frame bounds (cf. [2|), this procedure will result in 
a (compactly supported) shearlet frame for L^(i2) with the same frame bounds as 
before. 

We now apply this procedure to the family of compactly supported shearlet 
frames forL^(R^), which yielded (almost) optimally sparse approximations of func- 
tions in S'y (see ifTSl Thm. 1.3]). The main result of this paper then proves that the 
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resulting family of shearlet frames - now regarded as a system on [0, 1]^ with com- 
pact support in Q - again provides (almost) optimally sparse approximations now 
of elements from our model of cartoon-like images on bounded domains S"^ ^{Q) in 
the sense of 

\\f-fN\\2<C-N-'--{\ogNf asN-^oo. 

The precise statement is phrased in Theorem[T]in Section|3] 
This result is quite surprising in two ways: 

• Surprise 1. Regarding a log-factor as negligible - a customarily taken viewpoint 
-, the previous result shows that even for our much more sophisticated model of 
cartoon-like images on bounded domains the same optimal sparse approximation 
rate as for the simple model detailed in Subsection II. II can be achieved. This is 
even more surprising taking into account that our model contains point singular- 
ities at the corner points of the singularity curves. Naively, one would expect that 
these should worsen the approximation rate. However, observing that 'not too 
many' shearlets intersect these 'sparsely occurring' points unravels this mystery. 

• Surprise 2. Orthogonally projecting a shearlet system onto the considered boun- 
ded domain, thereby merely truncating it, seems an exceptionally crude approach 
to derive shearlets for a bounded domain. However, these 'modified' shearlet 
systems are indeed sufficient to achieve the optimal rate and no sophisticated 
adaptions are required, which is of significance for deriving fast algorithmic re- 
alizations. 



1.6 Main Contributions 

The main contributions of this paper are two-fold. Firstly, we introduce S^ ^[£1) 
as a suitable model for a function on a bounded domain with curvilinear singular- 
ities. Secondly, we show that the 'crude' approach towards a shearlet system on a 
bounded domain by simply orthogonally projecting still provides optimally sparse 
approximations of elements belonging to our model class S^ ^[Q.). 

We should mention that although not formally stated the idea of one piecewise 
C^ discontinuity curve in a model for functions on M^ as an extension of Defini- 
tion [T] is already lurking in 1 1 J . Also a brief sketch of proof of (almost) optimally 
sparse approximations of curvelets is contained therein. These ideas are however 
very different from ours in two aspects. First of all, our goal is a suitable model 
for functions on bounded domains exhibiting discontinuity curves and also treat- 
ing the boundary of the domain as a singularity curve. And secondly, in this paper 
we consider compactly supported shearlets - hence elements with superior spatial 
localization properties in contrast to the (band-limited) curvelets - which allows an 
elegant proof of the sparse approximation result in addition to a simplified treatment 
of the corner points. 
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1.7 Outline 

In Section|2] after recalling the definition of shearlet systems, we introduce shearlet 
systems on bounded domains, thereby focussing in particular on compactly sup- 
ported shearlet frames. The precise statement of our main result is presented in 
Section[3]together with a road map to its proof. The proof itself is then carried out 
in Section|4] Finally, in Section|5] we discuss our results and possible extensions of 
it. 



2 Compactly Supported Shearlets 

We first review the main notions and definitions related to shearlet theory, focussing 
in particular on compactly supported generators. For more details we would like to 
refer the interested reader to the survey paper 1,14]. Then we present our definition 
of shearlet systems on a bounded domain Q. S STAR^{v,L). 



2.1 Compactly Supported Shearlet Frames for L^{M?) 

Shearlets are scaled according to a parabolic scaling law encoded in the parabolic 
scaling matrices A2J orA2/, j G Z, and exhibit directionality by parameterizing slope 
encoded in the shear matrices Si^, k eZ, defined by 

,^2^' \ - fv/^ 

At, = ;/T or At, = 



and 



2J/y "^ "2. -^Q 2J 



^'~[o\ 



respectively. 

We next partition the frequency plane into four cones "^i - '^4. This allow the 
introduction of shearlet systems which treat different slopes equally in contrast to 
the shearlet group-based approach. We though wish to mention that historically the 
shearlet group-based approach was developed first due to very favorable theoretical 
properties and it still often serves as a system for developing novel analysis strate- 
gies (see, for instance, IfTTI ). 

The four cones "^i - "^4 are now defined by 

{i^u^2)&«':^i>lA^2/^i\<l}--i = l, 
{(<^,,<^2)6M':^i<-l,|&/<^i|<l}:i = 3, 
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and a centered rectangle 

^ = {(^i,&)gM2: 11(^1, ^2)||oo<l}. 
For an illustration, we refer to Figure |2a). 



'&2 



^3 




^1 



■Kt \^ 




(a) 



(b) 



Fig. 2 (a) The cones "^i - ^4 and the centered rectangle Sf, in frequency domain, (b) The tihng of 
the frequency domain induced by a cone-adapted shearlet system, where the (essential) support of 
the Fourier transform of one shearlet generator is exemplary high-lighted. 



The rectangle ^ corresponds to the low frequency content of a signal and is 
customarily represented by translations of some scaling function. Anisotropy comes 
into play when encoding the high frequency content of a signal which corresponds 
to the cones '^i\ - '^#4, where the cones '£1 and "^3 as well as '£2 and "^4 are treated 
separately as can be seen in the following 

Definition 4. For some sampling constant c > 0, the cone-adapted shearlet system 
SH{^,\l/,\(f;c) generated by a scaling function (j) G L?{M?) and shearlets Xj/jij/ £ 



) is defined by 



where 



SH{^,\l/,\lf;c) ^ 4>(0;c)Uf (v/;c)U'f'(v/;c), 
<P{<j);c) = {(j),„ = (j){--cm) : m e I?}, 



W{^if-c) = {v/,-.i,,„ - l^il^wiSkAy ■ -cm) : ;■ > 0, |^| < \Vl\m e I?}, 



and 



«F(va;c) = { v/,. ,,,„ = 23^/V(5[A2, • -cm) : ; > 0, |^| < \2'l^],m G I?}. 

The tiling of frequency domain induced by SH{(j),\i/,\(r;c) is illustrated in Fig- 
ure |2b). From this illustration, the anisotropic footprints of shearlets contained in 
f (v/;c) and ^{\[r;c) can clearly be seen. The corresponding anisotropic footprints 
of shearlets in spatial domain are of size 2^-' " times 2^-'. 
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The reader should keep in mind that although not indicated by the notation, the 
functions 0,,,, i/^.A,™- and Wj.k.m ^H depend on the sampling constant c. For the sake 
of brevity, we will often write y/i and xj/x, where A = {j,k,m) index scale, shear, 
and position. For later use, we further let Aj and Aj be the indexing sets of shearlets 
in f (\//;c) and 'f'(\J/;c) at scale j, respectively, i.e., 

f(V/;c) = {v/A:AeA,-,y = 0,l,...}and'f'(vA;c) = {v/;L:AeA,,y=0,l,...}. 

Finally, we define 

A=\jAj and A^[JAj. 

7=0 7=0 

The shearlet systems SH{<^,^r,^r\c) have already been very well studied with 
respect to their frame properties for L^(IR^), and we would like to refer to results 
in Q [l2l O. It should be mentioned that those results typically concern frame 
properties of f (i/a;c), which immediately imply frame properties of 'P{\fr;c) like- 
wise, whereas numerous frame properties for the low-frequency part 0(0;c) can 
be found in the wavelet literature. Combining those leads to frame properties of 
SH{(j),\i/,\i/;c). 

Recent results in fTOl establish frame properties specifically for the case of spa- 
tially compactly supported shearlet systems, i.e., shearlet systems with compactly 
supported generators (j), \j/, and iff which lead to a shearlet system consisting of 
compactly supported elements. These results give sufficient conditions for the so- 
called ttj conditions to be satisfied. As one class of examples with 'good' frame 
bounds, generating shearlets y/ and i/A were chosen to be separable, i.e., of the form 
V/'i(jci) ■ V'2(j^2) and '/ife) • V'2(-^i)i respectively, where y/i is a wavelet and \//2 a 
scaling function both associated with some carefully chosen (maximally flat) low 
pass filter. The separability has in addition the advantage to lead to fast accompany- 
ing algorithms. 

We wish to mention that there is a trade-off between compactly support of the 
shearlet generators, tightness of the associated frame, and separability of the shear- 
let generators. The known constructions of tight shearlet frames do not use separable 
generators, and these constructions can be shown to not be applicable to compactly 
supported generators. Tightness is difficult to obtain while allowing for compactly 
supported generators, but we can gain separability, hence fast algorithmic realiza- 
tions. On the other hand, when allowing non-compactly supported generators, tight- 
ness is possible, but separability seems to be out of reach, which makes fast algo- 
rithmic realizations very difficult. 



2. 2 Compactly Supported Shearlet Frames for l}{Q.) 

Let now Q. E STAR^{v,L) be a bounded domain as defined in Subsection [T3] The 
main idea to introduce a shearlet frame for L^{Q), preferably with compactly sup- 
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ported elements, is to start with a compactly supported shearlet frame forL'^(R^) and 
apply the orthogonal projection onto I?{Q,) to each element. To make this mathe- 
matically precise, we let Pq : L^{M?) -^ I?{Q.) denote the orthogonal projection 
onto L?-{Q.). 

Definition 5. Let Q e STAR'-{v,L). For some sampling constant c > 0, the cone- 
adapted shearlet system SHq ((^, y/, yf',c) for L (Q) generated by a scaling function 
(j) e L^(M^) and sliearlets \l/,yff & L^iR^) is defined by 

SHni^,Vf,T,c)^Pnm^;c)u^'{xi/;c)U^{xif-c)), 

where 0{^;c), f (v/;c), and 'P{\fr;c) are defined as in Definition^ 

As a direct corollary from well known results in frame theory (see Q), we obtain 
the following result, which clarifies frame properties for systems SHq (0 , i//, \J/; c) to 
the extent to which they are known for systems SH{(j), y/, ij/; c). In the sequel, we will 
usually regard SHq{(J),\i/, \ff;c) as a system defined on [0, 1]^ - in accordance with 
our model S^ ^{Q.) - by which we simply mean extension by zero. This system 
will be sometimes referred to as the extension of SHq{^,\i/,\j/;c) to [0,1]^. The 
following result also provides frame properties of these systems. 

Proposition 1. Let c> Q, let (^ ,\\f ,\\f ^ ^'(M^), and let Q. e STAR^{v,L) with pos- 
itive measure. Then the following conditions are equivalent. 

(i) The shearlet system SH{(p, Xj/, V/;c) is a frame for L (R. ) with frame bounds A 

and B. 
(ii) The shearlet system SHq {(j), y/, yf',c) is a frame for L (Q) with frame hounds A 

and B. 
(Hi) The extension of the shearlet system SHq{(J), y/, yf;c) to [0, 1] is a frame with 

frame bounds A and B for functions L" ( [0, 1 ] ) with compact support in Q. 



3 Optimal Sparsity of Shearlets on Bounded Domains 

We now have all ingredients to formally state the result already announced in Sub- 
section 11.51 which shows that even with the 'crude' construction of shearlets on 
bounded domains and the significantly more sophisticated model for cartoon-like 
images on bounded domains we still obtain (almost) optimally sparse approxima- 
tions. 

Theorem 1. Let c > 0, and let (p,y/,y/ E L~{M. ) be compactly supported. Suppose 
that, in addition, for all ^ = (^i , ^2) G IK^'. the shearlet yr satisfies 



(i) 
(ii) 



1^(^)1 <Ci-min(l,|^i|«).mina,|<^ir^)-min(l,|fe|-r),W 
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where a > 5, y > A, h E L (R), and C\ is a constant, and suppose that the shearlet 
iff satisfies (i) and (ii) with the roles of^\ and ^2 reversed. Further, letv>Q,LG Z+ 
and Q € STAR {v,L), and suppose that SHq((J), y/, \i/;c) forms a frame for L (Q). 
Then, the extension of the shearlet frame SHQ{(p,y/,\j/;c) to [0,1] (cf Sub- 
section \2.2]l provides (almost) optimally sparse approximations of functions f G 
(oy i^{Q) in the sense that there exists some C > such that 

\\f-fN\\l<C-N-^-i\ogN)^ asN^c., 

where f^ is the nonlinear N-term approximation obtained by choosing the N largest 
shearlet coefficients off. 



3.1 Architecture of the Proof of Theorem\l\ 

Before delving into the proof in the following section, we present some preparation 
before as well as describe the architecture of the proof for clarity purposes. 

Let now SHQ[^,y/,\jf;c) satisfy the hypotheses in Theorem [T] and let / G 
S'y ^{Q.). We first observe that, without loss of generality, we might assume the 
scaling index j to be sufficiently large, since / as well as all frame elements in the 
shearlet frame SHQ{(p,\j/,y/;c) are compactly supported in spatial domain, hence 
a finite number does not contribute to the asymptotic estimate we aim for In par- 
ticular, this means that we do not need to take frame elements from 0(0;c) into 
account. Also, we are allowed to restrict our analysis to shearlets y/j,k,m, since the 
frame elements yfj,k.m can be handled in a similar way. 

We further observe that we can drive the analysis for the frame SH{(p,\i/, yf;c) 
and for the domain [0, 1]^ instead, since, by hypothesis, Q is contained in the interior 
of [0, 1]^, we treat the boundary of i2 as a singularity curve in [0, 1]^, and the frame 
properties are equal as shown in Proposition [T] In this viewpoint, the function to be 
sparsely approximated vanishes on [0, 1]^ \I2. 

Our main concern will now be to derive appropriate estimates for the shearlet co- 
efficients {{f,\i/i):XGA}off. Letting 1 (/) |„ denote the «th largest shearlet coef- 
ficient (/, yf^ ) in absolute value and exploring the frame property of ^//(^ , y/, iff; c), 
we conclude that 

ll/-/iv||?<^E|0(/)l^, 

^ n>N 

for any positive integer A^, where A denotes the lower frame bound of the shearlet 
frame SH{^ , i//, iff; c) . Thus, for the proof of Theoremfl] it suffices to show that 

Y,\0{f)\l<C-N-^-ilogN)' asN^^. (3) 

n>N 

To derive the anticipated estimate in Q, for any shearlet yf^, we will study two 
separate cases: 
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• Case 1. The compact support of the shearlet ^^x does not intersect the boundary 
of the setB(or5f2), i.e., supp(v/;t) n (5BU5i2) =0. 

• Case 2. The compact support of the shearlet i//;l does intersect the boundary of 
the set B (or (9 f2), i.e., supp(v/A) n (<9BU(9i2) ^0. 

Notice that this exact distinction is only possible due to the spatial compact support 
of all shearlets in the shearlet frame. 

Case 2 will then throughout the proof be further subdivided into the situations 
- which we now do not state precisely, but just give the reader the intuition behind 
them: 

• Case 2a. The support of the shearlet does intersect only one C^ curve in dB U dO. . 

• Case 2b. The support of the shearlet does intersect at least two C^ curves in 
dB\JdQ.. 

- Case 2b- 1. The support of the shearlet does intersect dBUdQ in a corner 
point. 

- Case 2b-2. The support of the shearlet does intersect two C^ curves in dB U 
dQ simultaneously, but does not intersect a corner point. 



4 Proof of Theorem [T] 

In this section, we present the proof of Theorem [H following the road map out- 
lined in Subsection 13. II We wish to mention that Case 1 and Case 2a are similar 
to the proof of (almost) optimally sparse approximations of the class <o'^(v) using 
compactly supported shearlet frames in ifTSJI . However, Case 2b differs significantly 
from it, since it, in particular, requires a careful handUng of the corner points of dB 
and dQ . 

In the sequel - since we are concerned with an asymptotic estimate - for sim- 
pUcity we will often simply use C as a constant although it might differ for each 
estimate. Also all the results in the sequel are independent on the sampling constant 
c > 0, wherefore we now fix it once and for all. 



4.1 Case 1: The Smooth Part 

We start with Case 1, hence the smooth part. Without loss of generality, we can 
consider some g e C^([0, 1]^) and estimate its shearlet coefficients. The following 
proposition, which is taken from |15| , implies the rate for optimal sparsity. Notice 
that the hypothesis on \// of the following result is implied by condition (i) in Theo- 
rem [T] 

Proposition 2 dHU). Let g G C^([0, 1]^), and let Xj/ E L^(R^) be compactly sup- 
ported and satisfy 
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|r(^)|<Ci -111111(1, |^ir).min(l,|^i|-^).min(l, 1^21-^) /or flW^ = (^i, ^2) eM2, 

where 7 > 3, a > 7+ 2, and C\ is a constant. Then, there exists some C > such 
that 

Y^\e{g)\l<C-N'^ asN ^00. 

n>N 

This settles Theorem[T]for this case. 

4.2 Case 2: The Non-Smooth Part 

Next, we turn our attention to the non-smooth part, and aim to estimate the shear- 
let coefficients of those shearlets whose spatial support Intersects the discontinuity 
curve dB or the boundary of the domain Q. One of the main means of the proof will 
be the partitioning of the unit cube [0, l]' into dyadic cubes, picking those which 
contain such an intersection, and estimating the associated shearlet coefficients. For 
this, we first need to introduce the necessary notational concepts. 

For any scale j>Q and any grid point p E 1?, we let Qjp denote the dyadic cube 
defined by 

Further, let Qj be the collection of those dyadic cubes Qjp which intersect dBiJdQ, 
i.e., 

Qj = {Qj.p ■■ Qj,pn{dBUdQ)y^(d,p El}}. 

Of Interest to us is also the set of shearlet indices, which are associated with shearlets 
Intersecting the discontinuity curve inside some Qjp E Qj, hence, for j >Q and 
p El? with Qjp E Qj, we will consider the index set 

Aj^p = {Xe Aj : supp(v/;i) n Qj^p n{dBLidn)y^ 0}. 

Finally, for j > 0, p E 1?, and < e < 1, we define Ay,p(e) to be the index set 
of shearlets i//;^' ■^ G ^j,/" ^'^'-^ that the magnitude of the corresponding shearlet 
coefficient (/, y^x ) is larger than e and the support of i//^ Intersects Qj^p at the jth 
scale, i.e., 

A>(e) = {^eA,,p:|(/,v/A)|>e}, 

and we define A(e) to be the index set for shearlets so that \{f,'^fx)\ > £ across all 
scales j, i.e., 

A(e) = UA,-,,(e). 
j.p 

The expert reader will have noticed that in contrast to the proofs in H] and ||9l, which 
also split the domain into smaller scale boxes, we do not apply a weight function to 
obtain a smooth partition of unity. In our case, this is not necessary due to the spatial 
compact support of the frame elements. Finally we set 
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Sj-P^ U supp(v/a), 

which is contained in a cubic window of size C ■ 2^^'^ by C ■ 2^-'' ^, hence is of 
asymptotically the same size as Qj^p- As mentioned earlier, we may assume that j 
is sufficiently large so that it is sufficient to consider the following two cases: 

• Case 2a. There is only one edge curve Fi C dB (or dQ) which can be parame- 
terized by x\ ~ E{x2) (orx2 — E{xi)) with E gC^ inside 5/, p. For any A G Ajp, 
there exists somei = {xi^xi) G Qj,pr\supp{\i/^) CiFi. 

• Case 2b. There are two edge curves ri,l2 C dB (or dQ) which can be pa- 
rameterized by x\ ~ E{x2) (or X2 ~ E{x\)) with E G C^ inside Sj,p. For any 
A G Aj^p, there exist two distinct points x = (ii,i2) and y — {yi,y2) such that 

X e Qj,p n supp( v/A ) n Fi and y G Qj,p n supp( i/A ) n r2. 

In the sequel, we only consider the edge curve dB to analyze shearlet coefficients 
associated with the non-smooth part, since the boundary of the domain Q can be 
handled in a similar way; see also our elaboration on the fact that WLOG we can 
consider the approximation on [0, 1]^ rather than Q in Subsection l3.1l 



4.2.1 Case 2a: The Non-Smooth Part 

This part was akeady studied in ITSl . where an (almost) optimally sparse approx- 
imation rate by the class of compactly supported shearlet frames 5//(0, i//, V/;c) 
under consideration was proven, and we refer to |;15 1 for the precise argumentation. 
For intuition purposes as well as for later usage, we though state the key estimate, 
which implies (almost) optimally sparse approximation for Case 2a: 

Proposition 3 (||15|). Let Hf G L (R ) be compactly supported and satisfy the con- 
ditions (i), (ii) in Theorem\l\and assume that, for any A G Ajp, there exists some 
X = (ii , X2 ) G Qj,pf^ supp ( VX ) 1^ ^B. Let s be the slop^ of the tangent to the edge 
cun'e dB at (ii,i2), i-^-, 

• s = E'[x2), if dB is parameterized by x\ = E{x2) with E G C in Sj^p, 

• s = {E' {x\))^ , if dB is parameterized by X2 — E{xi) with E (^ C in Sj,p, and 

• s ^ °o, if dB is parameterized by X2 = E{xi ) with E'{xi ) = and E E C in Sj^p. 

Then, there exists some C > such that 

|(/,V/a)|<C-2-^\ if\s\>^or\s\^o.^ (4) 

and 



' Notice that here we regard the slope of the tangent to a curve {E{x2),X2), i.e., we consider i of 
a curve xi = SX2 + b, say. For analyzing shearlets V/y/.m, the roles of Xi and X2 would need to be 
reversed. 



Shearlets on Bounded Domains 



15 



\{f,Wi)\<C- 



2-^J 



if\s\<3. 



(5) 



|/t+2i/2i|3' 

Similar estimates with dB substituted by dQ hold if, for any X E A^ p, there 
exists some x = (^1,^2) G Qj,p H supp(v/x) n dQ. 



4.2.2 Case 2b: The Non-Smooth Part 

Letting e > 0, our goal will now be to first estimate |Aj^p(e)| and, based on this, 
derive an estimate for |A(e)|. WLOG we might assume ||)/||i < 1, which implies 



l(/,rA>l<2- 



5^ 



Hence, for estimating \Aj p(e)|, it is sufficient to restrict our attention to scales j < 
|log2(e-i). 

As already announced before, we now split Case 2b into the following two sub- 
cases: 

• Case 2b- 1. The shearlet i//^ intersects a corner point, in which two C^ curves Fi 
and I2, say, meet (see Figure[3](a)). 

• Case 2b-2. The shearlet VX intersects two edge curves Fi and I2, say, simultane- 
ously, but it does not intersect a corner point (see Figure |3](b)). 





(a) 



(b) 



Fig. 3 (a) A shearlet i/z^ intersecting a comer point where two edge curves Fi and F meet. T\ and 
T2 are tangents to the edge curves F and F in this comer point, (b) A shearlet i//^ intersecting two 
edge curves F and F which are a part of the boundary of sets Bq and B\.T\ and T2 are tangents to 
the edge curves F and F in points contained in the support of Xfr^ ■ 



Case 2b- 1. We first consider Case 2b- 1. In this case, by a counting argument, it 
follows that 

|A,,p(£)|<C-2^'/2. 

Since there are only finitely many corner points with its number not depending on 
scale j > 0, we have 
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|A(e)|<C- £ 2>/2<C-e-5. 

The value e > can be written as a function of the total number A^ of coefficients, 
which yields e(A^) < C-N^T-. This implies that 

l^\e[f)\l<c-N-\ 

n>N 

and the optimal sparse approximation rate is proven for Case 2b- 1. 

Case 2b-2. Next, we consider Case 2b-2. In this case, WLOG, we might assume 
that, for any A G A^ p, there exist two distinct points x = (ii ,X2),y — {yi,y2) such 
thatie Qj,pr\supp{\i/x)r\ri andy E Qj_pnsupp{\i/i)r\r2, and the two edge curves 
I] and I2 are parameterized by xj = E{x2) (or X2 = E{xi)) with E eC^ inside Sj^p. 
We can then write the function / G (#^^(i2) as 



foXBo +/1ZB1 = ifo-fi)XBQ+fi on Sj.p, 

where /o,/i E C^([0, 1]^) and Bq,Bi are two disjoint subsets of [0, 1]^ (see Figure 
O. By Proposition |2l the rate for optimal sparse approximation is achieved for the 
smooth part /i . Thus, it is sufficient to consider / — goXBg with go ~ fo — fi G 

C2([0,l]2). 

Assume now that the tangents to the edge curves Fi and I2 at the points x and y 
are given by the equations 

Ti : xi+ai=si{x2 + bi) and T2 : xi+a2 = S2{x2+b2), 

respectively, i.e., ^i and ^2 are the slopes of the tangents to the edge curves Fj and 
^2 at i and y, respectively. If the curve f;, / = 1,2, is parameterized by X2 = E{xi) 
with E'{xi) = 0, we let s, = 00 so that the tangent is given by X2 = —bi in this case. 
Next, for fixed scale j and shear index k, let N^j i^{Qj,p) denote the number of 
shearlets y/^ intersecting rj in Qj,p, i.e., 

NjAQj.p) = |{A = U,k,m) : Qj„ nsuppiy/^) nFi ^ 0}|, 

let A^?^(2j p) denote the number of shearlets Xj/i intersecting I2 in Qj.p^ i-e-- 

Nlk{Qj,p) = |{A = U,k,m) : Qj„nsupp{xif^)nr2 + 0}|, 

and let Nj^Qj^p) denote the number of shearlets i//;l intersecting Fj and I2 in Qi,p^ 
i.e., 

Nj,k{Qj,p) - |{A = (;,^,m) :ey,pnsupp(v/A)nri ^0 and e,-,pnsupp(v/;i)nr2 ^0}|. 

Then 

A^;,^(e,,/.) < mm{N]AQj,p),Nl,{Qj,p)). (6) 
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By a counting argument, there exists some C > such that 

Nl,{Qj,p)<C-V'^ for/- 1,2, (7) 

and the form of supp(v/)L) impHes 

N)^k{Qj.p) < C ■ {\2J/-st+k\ + 1) for / = 1,2. (8) 

We now subdivide into three subcases, namely, l^i | , 1^2 1 < 2, and \si\ < 2, |i2 1 > 2 
(or vice versa), and \si\, \s2\ > 2, and show in each case the (almost) optimal sparse 
approximation rate claimed in Theorem[T] This then finishes the proof. 

Subcase \si\, \s2\ < 2. In this case, (|6]l and ^ yield 

NjAQj.p) < C-min(|2^'/2ii+^| + l,\2J/\+k\ + 1). 

We first show independence on the values of si and S2 within the interval [—2,2]. 
For this, let s and s' be the slopes of the tangents to the edge curve Fi (or I2) at 
t S Qj.p n supp(i//x) and f' G Qj.,, D supp{\j/ii), respectively, with s G [—2,2]. Since 
Fi (or I2) is C^, we have 1^ — .?'| < C • 2^J'^, and hence 

\2J/^s'+k\<C-{\2J/^s+k\ + l). 

This implies that the estimate for Nj,k{Qj.p) asymptotically remains the same, in- 
dependent of the values of si and S2. Further, we may assume s' G [—3,3] for 
s G [—2,2], since a scaling index j can be chosen such that [>?- i'| is sufficiently 
small. Therefore, one can apply inequality Q from Proposition [3] for both points t 
and t'. In fact, it can be easily checked that one can use ^ with the slope s instead 
of s' for the point t' (or vice versa); this replacement will not change the asymptotic 
estimates which we will derive. Thus, we might from now on use universal values 
for the slopes si and S2 at each point in Qj,p- 
Now using (|5]l from Proposition[3] we have 

|(/,V/,M.)I <C-max(-|£-,-Ai_). (9) 



2J/hi+k\^'\2J/\^2+k\^ 



Since X ^ ,,3 > £ implies 



\2J/^Si+k\ < e-^2-^J for / = 1,2, 
the estimate (|9|l yields 
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|A;,p(e)| <C- Y. min(|2^'/2^i+^| + l,|2-''/2i2+^| + l) 

keK]iE]UKJ(E) 
>=lkeK'j{E) 

<c-{£-h-lJ + i)^, (10) 

where 

K'j{e) ^{keZ: \2'/^si+k\ < C-e^h^^J} for / = 1,2. 

By the hypothesis for Case 2b-2, we have \Qj\ < C, where the constant C is inde- 
pendent of scale j > 0. Therefore, continuing ( fTOl i. 

llog2(£"') 

\A{e)\<C- £ |A,-,,(e)|<C.e-i 

7=0 

This allows us to write £ > as a function of the total number of coefficients A^, 
which gives 

e(A^) <C-N-i. 

Thus 

Y,\eif)\l<c-N-\ (11) 

n>N 

which is the rate we seeked. 

Subcase |si | < 2 and 1^2 1 > 2 or vice versa. In this case, ©-([SI) yield 

Nj,k{Qj.p)<C-min{\2J/^si+k\ + l,2''^)- 

Again utilizing the fact that the edge curves are C^, and using similar arguments as 
in the first subcase, WLOG we can conclude that the slopes si,S2 at each point in 
Qj,p are greater than j. 



Now, exploiting inequalities (|4]i and (|5]l from Proposition[3] we have 

_ 3 . 

\{f,WjX,n)\<C-m.^( J'' 2-'^^). (12) 



2J/hi+k\^' 
\2^^^si+k\ <e'h-^J, 



Since |,,/1.,''_^,|3 > £ implies 



9 . 



and 2 3^ > £ implies 
it follows from (fT2T i. that 



;<^log2(£-'), 
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3log2(£-') |log2(<:-') 

.7=0 k&KJ(£) j=0 

<C-( £ (e-h-^'/4 + i)2+ ^ 2^'/2 
;=o 7=0 

The value e > can now be written as a function of the total number of coefficients 
A^, which gives 

e{N) <C-N-i. 

Thus, we derive again the seeked rate 



n>N 

Subcase \si\ > 2 and \s2\ > 2. In this case, (|6]l and dTji yield 
Following similar arguments as before, we again derive the seeked rate (fTTl i. 

5 Discussion 

A variety of applications are concerned with efficient encoding of 2D functions 
defined on non-rectangular domains exhibiting curvilinear discontinuities, such as, 
e.g., a typical solution of a transport dominated partial differential equation. As an 
answer to this problem, our main result, Theorem[T] shows that compactly supported 
shearlets satisfying some weak decay and smoothness conditions, when orthogo- 
nally projected onto a given domain bounded by a piecewise C^ curve, provide 
(almost) optimally sparse approximations of functions which are C^ apart from a 
piecewise C' discontinuity curve. In this model the boundary curve is treated as a 
discontinuity curve. 

Analyzing the proof of Theoremfl] it becomes evident that the presented optimal 
sparse approximation result for functions in S'yj^{Q) generalizes to an even more 
encompassing model, which does contain multiple piecewise C" possibly intersect- 
ing discontinuity curves separating C^ regions in the bounded domain Q . 

In some applications it is though of importance to avoid discontinuities at the 
boundary of the domain. Tackling this question requires further studies to carefully 
design shearlets near the boundary, and this will be one of our objective for the 
future. 
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